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The impact of quenched disorder on deterministic diffusion in chaotic dynamical systems is studied. 
As a simple example, we consider piecewise linear maps on the line. In computer simulations we find 
a complicated scenario of multiple suppression and enhancement of normal diffusion under variation 
of the perturbation strength. These results are explained by a theoretical argument showing that 
the oscillations emerge as a direct consequence of the unperturbed diffusion coefficient, which is 
known to be a fractal function of a control parameter. 
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Recently there has been growing interest in the field 
of disordered dynamical systems thus trying to bring 
together two at first view very different directions of 
research Disordered lattices exhibiting quenched 

(static) randomness are considered as a traditional prob- 
lem of statistical physics. Hence, in this type of models 
normal and anomalous diffusion are studied by proba- 
bilistic methods being developed in the framework of the 
theory of stochastic processes On the other hand, 

diffusion can also be generated from microscopic chaos 
in nonlinear equations of motion. This is well-known as 
the phenomenon of chaotic diffusion in deterministic dy- 
namical systems Here, methods of dynamical sys- 
tems theory can be applied for computing deterministic 
transport coefficients ||^-p^. Disordered dynamical sys- 
tems represent an interesting combination of both types 
of models and provide an ideal opportunity to bring these 
two theoretical approaches together. To our knowledge, 
only very few cases of such models have been studied so 
far. Examples include random Lorentz gases for which 
Lyapunov exponents have been calculated by means of 
kinetic theory and by computer simulations [0|, numer- 
ical studies of diffusion on disordered rough surfaces p^ ] 
and in disordered deterministic ratchets jist , as well as 
numerical and analytical studies of chaotic maps on the 
line with quenched disorder P,p^. 

In this work we will focus on the most simple case in 
the latter class of models, which are piecewise linear maps 
defined on the unit interval and acting on the real line 
by periodic continuation. In case of mixing dynamics, 
the unperturbed maps exhibit normal diffusion |f7|,p^-^ . 
As first shown in Refs. [^, and as further analyzed in 



Refs. |14|, quenched disorder can change the determinis- 
tic diffusive dynamics in these maps profoundly: adding 
static randomness in form of a local bias with globally 
vanishing drift leads to dynamical localization of trajec- 
tories in a complicated potential landscape. On disor- 
dered lattices, this effect is well-known as the Golosov 
phenomenon [p|p9|] thus reappearing in the framework of 
deterministic dynamics. Here we will consider a second 
fundamental type of static randomness which is multi- 
plicative and preserves the local symmetry of the model. 
Consequently, it does not generate subdiffusive behavior 
ipof . On this basis we study the dependence of the normal 



diffusion coefficient on two control parameters, which are 
the strength of the unperturbed diffusion coefficient as 
well as the perturbation strength [pl| . Under variation 
of these two control parameters we find a complicated 
scenario for the diffusion coefficient yielding suppression 
and enhancement of the strength of diffusion as general 
features. These findings will be explained in terms of a 
simple theoretical approximation, which is derived from 
an exact diffusion coefficient formula as obtained within 
the theory of stochastic processes. 

We first define the unperturbed model by the equation 
of motion 
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FIG. 1. Diffusion coefficient D{a,da) for the piecewise lin- 
ear map shown in the figure. The slope a is perturbed by 
static disorder of maximum strength da as defined in Eq. ^ . 
The bold black line depicts numerically exact results for the 
unperturbed diffusion coefficient at da = 0. Computer simu- 
lation results for da > are marked with symbols, the corre- 
sponding lines are obtained from the theoretical approxima- 
tion Eq. .(P) The parameter values are: da — 0.1 (circles), 
da — 0.4 (squares), da = 1.0 (diamonds). 

where a £ M is a control parameter and Xn is the posi- 
tion of a point particle at discrete time n. Ma{x) is con- 
tinued periodically beyond the interval [—1/2,1/2) onto 
the real fine by a lift of degree one, Ma{x+1) = Ma(a:)-I-1. 
We assume that Ma{x) is anti-symmetric with respect to 
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a; = 0, Ma{x) = —Ma{—x). The map we study as an 
example is defined by Ma{x) — ax, where the uniform 
slope a serves as a control parameter. The Lyapunov 
exponent of this map is given by A(a) = Ina imply- 
ing that for a > 1 the dynamics is chaotic. We now 
modify this system by adding a random variable Aa(i), 
J G Z, to the slope on each interval [i — 1/2, i + 1/2) 

yielding Ma\x) = (a + ha{i))x. We assume that the 
random variables Aa(i) are independent and identically 
distributed according to a distribution Xda(Aa), where 
da is again a control parameter. In the following we will 
consider two different types of such distributions, namely 
random variables distributed uniformly over an interval 
of size [—da, da], 

Xda(Aa) = {Q{da + Aa)e(da - Aa))/(2da) , (2) 

and dichotomous or ^-distributed random variables, 

X<ia(Aa) ((5(da-}- Aa) +(5(da- Aa))/2 . (3) 

Since |Aa| < da, we denote da as the perturbation 
strength. As an example, we sketch in Fig. |l| the map 
resulting from the disorder of Eq. (H) as applied to the 
slope a = 3. In the absence of any bias, the diffusion co- 
efficient is defined as D{a,da) = lim„^oo < a:^ > /(2n), 
where the brackets denote an ensemble average over mov- 
ing particles for a given configuration of disorder. An 
additional disorder average is not necessary because of 
self-averaging. Note that for locally symmetric quenched 
disorder and (a — da) > 2 there is no physical mechanism 
leading to infinitely high reflecting barriers as they are 
responsible for Golosov localization ||^. Thus diffusion 
must be normal, as is confirmed by computer simula- 
tions. Hence, the central question is what happens to the 
parameter-dependent diffusion coefficient D(a, da) under 
variation of the two control parameters a and da. 

For the unperturbed case da = the diffusion coef- 
ficient has been computed numerically exactly in Refs. 
[^,0. There it has been shown that D{a,0) is a frac- 
tal funtion of the slope a as a control parameter. This 
function is depicted in Fig. 0, as well as results from 
computer simulations for different values of the perturba- 
tion strength da [Q. As expected, the fractal structure 
smoothes out by increasing da. However, it is remark- 
able that even for large perturbation strength da oscil- 
lations are still visible as a function of a indicating that 
the original irregularities are very robust against uniform 
quenched disorder. 

Before we proceed to more detailed simulation results 
we briefiy repeat what is known for diffusion in lat- 
tice models with random barriers p|-p|,p3[. In the most 
simple version, the quenched disorder is defined on a 
one-dimensional periodic lattice with transition rates be- 
tween neighboring sites i and i + 1 having the symmetry 
= = Tfc for a given random distribution of 

Ffc . In this situation an exact expression for the diffusion 
coefficient has been derived reading P,p3[ 



D = {l/T}-H^ , (4) 

with the brackets defining the disorder average {l/F} = 
l/iV^^gl/Ffe at chain length N, and for a distance 
I between sites. The double-inverse accounts for the 
physical significance of the fact that the highest barri- 
ers dominate diffusion in d = 1. In other words, the 
existence of rates with F^ ^ naturally leads to a van- 
ishing diffusion coefficient D 0. This scenario is trans- 
lated to the map under consideration as follows: Eq. (|l|) 
can be understood as a time-discrete Langevin equation, 
Xn+l — Xji - dV/dx{xn) In case of quenched disor- 

der, and under proper integration of Ma\x), the result- 
ing potential V{x) corresponds to the one of a random 
barrier model in which the perturbation strength da de- 
termines the highest barriers. For the disordered map one 
may thus naively expect suppression of diffusion reflected 
in I? (a, da) being a monotoneously decreasing function of 
da. 
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FIG. 2. Diffusion coefficient D{a,da) as a function of the 
perturbation strength da at slope a = 6: (a) disorder distri- 
bution Eq. (b) disorder distribution Eq. (^. The circles 
represent results from computer simulations, the lines are ob- 
tained from the theoretical approximation Eq. 

To check this hypothesis, we choose fixed values of a 
corresponding to the two extreme situations of sitting at 
a local maximum or minimum, respectively, of the un- 
perturbed D{a,0) in Fig. 0. We first focus on the lo- 
cal minimum at a = 6 for da < 0.5 with uniform (Fig. 
§ (a)) as well as dichotomous (Fig. ^ (b)) disorder. In 
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sharp contrast to the stochastic argument outhned above, 
in both cases we observe enhancement of diffusion as a 
function of da. Moreover, this enhancement does not ap- 
pear in form of a simple functional dependence on da: 
In (a), smoothed-out oscillations are visible on smaller 
scales, whereas in (b) the resulting function is clearly 
non-monotoneous and wildly fluctuating exhibiting mul- 
tiple suppression and enhancement in different parameter 
regions. What happens on larger scales of da is depicted 
in Fig. 1^, where we only show results for dichotomous dis- 
order. As can be seen in Fig. ^ (b), choosing a at a local 
maximum of D(a, 0) leads to suppression of diffusion for 
da < 1.0, whereas a local minimum generates enhance- 
ment in the same parameter region of da. However, in 
both cases the diffusion coefficient decreases on a larger 
scale thus recovering qualitative agreement with the ex- 
pectation from stochastic theory. Indeed, for {a— da) 2 
barriers are formed which a particle cannot cross anymore 
implying the existence of localization [e3 . 



cient on the unperturbed lattice we rewrite Eq. as 
D — {l/£)(rfc, In case of our map, the transi- 

tion rates and the distance between sites are both some- 
what combined in the action of the slope a as a control 
parameter. Therefore, the unperturbed diffusion coeffi- 
cient is correctly rewritten by replacing D{Tk, I) = D{a). 
Up to this point we performed purely formal manipu- 
lations. The key question is now what function shall 
be used for the parameter-dependent diffusion coefficient 
D{a) in case of deterministic dynamics. Here we pro- 
pose to identify the function D{a) with the exact, unper- 
turbed deterministic diffusion coefficient previously de- 
fined as D(a,0). Providing this information, the exact 
formula Eq. for stochastic dynamics becomes an ap- 
proximation which can straightforwardly be applied to 
deterministic dynamics in disordered systems. If the dis- 
order distributions Xdai^o.) are bounded by the pertur- 
bation strength da, and if we go into the continuum limit 
for the random variable, our final result reads 




Dapp{a, da) = 



FIG. 3. Diffusion coefficient D{a, da) as a function of the 
perturbation strength da for dichotomous disorder Eq. (^) at 
two different slopes a: (a) a = 6, (b) a — 7. The circles repre- 
sent results from computer simulations, the lines are obtained 
from the theoretical approximation Eq. (^. 

To theoretically explain these simulation results, we 
modify Eq. (Q) in a straightforward way such that it can 
be applied to our disordered deterministic map under 
consideration. We first note that for uniform transition 
rates Ffc = const, it is D{Tk,l) = TkP. Using this fa- 
miliar expression for the random walk diffusion coeffi- 



" j^A ^ Xda(Aa) 

diAa) — T — 7 

da Dia + Aa) 



(5) 



This expression represents the central formula of our let- 
ter. The results obtained from it are depicted in Figs. 
|l| to H in form of lines. For small enough da the agree- 
ment between theory and simulations is excellent, both 
for dichotomous as well as for uniform disorder. Surpris- 
ingly, for dichotomous disorder the theory even correctly 
predicts the oscillations for larger da, although there are 
clear quantitative deviations. 

We now show that this formula provides a simple phys- 
ical explanation for the complicated dependence of the 
diffusion coefficient on the perturbation strength. In case 
of da 0, Taylor expansion leads to 



da 



Dapp{a, da) = / d{Aa) Xda{Aa)D{a + Aa) 



(6) 



da 



We remark that, alternatively, Eq. (g) can be proven for 
small enough da starting from the precise definition of 
the diffusion coefficient, and without employing Eq. (|^) 
|p5| . In this limit the perturbed diffusion coefficient re- 
duces to an average of the exact diffusion coefficient over 
the neighborhood in the parameter interval [a— da, a-\-da\ 
weighted by the respective disorder distribution XdaiAa). 
Consequently, if a is chosen at a local minimum the result 
must be enhancement of diffusion by increasing the per- 
turbation strength, and suppression at a local maximum, 
respectively [Q. On these grounds it is also clear that 
the fractal parameter dependence of D{a, 0) must reap- 
pear in the perturbed diffusion coefficient thus leading to 
multiple suppression and enhancement on fine scales. In 
case of dichotomous noise, and if Eq. holds, D{a, da) 
must be a superposition of two fractals resulting in a new 
fractal. In case of uniform noise the fractality smoothes 
out, but is still responsible for the oscillations on a fine 
scale. 
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We conclude with a few remarks: (a) It would be im- 
portant to have a proof of Eq. (H) for dynamical sys- 
tems, as well as to obtain higher order corrections for 
explaining the deviations between simulation and theory 
as visible in Fig. ^. (b) Our results might be important 
to understand diffusion on a stepped surface with a dis- 
ordered arran gem ent of Ehrlich-Schwoebel barriers: As 
shown in Ref . , this problem can be modeled by a po- 
tential relief consisting of a combination of random traps 
and random barriers. Our map provides a deterministic 
generalization of such a model and particularly enables to 
study the impact of long-range memory effects on surface 
diffusion, (c) One could think of applying our approach 
to systems such as the ones studied in Refs. |l^Jl^, or to 
the periodic Lorentz gas [§-0, which is a model being 
close to experiments on antidot lattices Knowing 
the density-dependent diffusion coefRcient in the unper- 
turbed case leads us to predicting local and global 
suppression of the diffusion coefficient in this model in 
case of static density fluctuations. 

We wish to thank Prof. K.W. Kehr for pointing out to 
us the existence of Eq. (|^), unfortunately, he did not live 
to see of how much use his hint was for the present work. 
The author is also grateful to G. Radons, H. van Bei- 
jeren, J.R. Dorfman, and T. Tel for helpful discussions. 
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